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developed  using  differential  phase  data  from  the  ATS-6  geostationary  satellite. 
The  40-  and  3 60 -MHz  coherent  signals  were  propagated  through  the  iono- 
sphere and  received  on  the  150-ft  radio  telescope  at  the  Sagamore  Hill  Radio 
Observatory,  v 

Thf  rnfr*^if7  of  the  computerized  data  was  based  on  13  consecutive  days 
scintillation  in  November  1974.  It  shows  that  amplitude  and  phase 
. scintillation  occurred  simultaneously  only  53  percent  of  the  time,  indicating 
/ that  application  of  the  equivalent  thin  diffracting  screen  theory  is  invalid  for 
^ describing  the  data  for  the  remaining  47  percent  of  the  time.  Statistical 
comparisons  suggest  that  daytime  amplitude  scintillations  were  apparently 
caused  by  ionospheric  focusing  mechanisms  since  phase  scintillation  proved 
to  be  almost  exclusively  a nighttime  phenomena. 

The  differential  phase  data  from  the  same  source  is  analyzed  in  a manner 
that  reveals  many  of  the  ionospheric  interactions  that  affect  space  communica- 
tions. The  interactions  are  further  reduced  into  the  following  parameters 
that  describe  ionospheric  conditions  that  relate  to  rf  propagation: 

1.  Total  electron  content  (TEC) 

2.  Changes  in  electron  density,  assuming  height 

3.  Equivalent  slab  thickness,  assuming  height 

4.  A continuous  measurement  of  phase  velocity 

5.  A continuous  measurement  of  rate  of  change  of  TEC 

6.  A continuous  measurement  of  rate  of  change  of  refractive  index 

7.  A continuous  measurement  of  angle  of  arrival. 

The  influence  that  the  traveling  ionospheric  disturbances  (TID)  have  on 
the  phase  of  the  40-MHz  signal  receives  special  emphasis,  and  it  is  shown 
how  the  presence  of  TID  elds  in  the  analysis. 
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ATS-6  40-  and  360-MHz 
Differential  Phase  Measurements 


1.  INTRODUCTION 

Detection  and  processing  of  VHF  and  UHF  Satellite  cw  signals  with  power 
levels  less  than  -135  dbm  at  the  antenna  present  difficulties  when  attempting  to 
translate  the  processed  data  into  meaningful  measurements  of  ionospheric 
physical  interactions.  The  principal  problem  is  interference  from  extraneous 
signals,  that  is,  man-made  and  natural  rf  noise.  For  several  years  tracking 
filters  have  been  employed  to  alleviate  the  problem  by  narrowing  the  receiver 
bandwidth  to  a few  cycles  per  second.  The  filters  are  required  to  frequency  track 
the  rf  signal,  because  generally  the  satellite  signal  and  the  local  oscillator  are 
not  sufficiently  stable  to  allow  use  of  a narrowband  fixed  filter. 

Experience  has  shown  that  geostationary  satellite  UHF  signals  drift  up  to 
200  Hz  and  orbiting  satellite  signals,  due  to  the  added  doppler  frequencies,  shift 
many  tirpes  that  amount.  Local  oscillators,  however,  can  be  made  extremely 
stable  using  synthesizer  techniques. 

The  complexity  of  the  tracking  filter  with  the  automatic  search  circuitry 
required  to  lock  the  filter  onto  the  signal  during  acquisition,  or  when  it  reappears 
after  a fade,  is  an  unattractive  feature  of  this  process. 


(Received  for  publication  1 June  1976) 
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Experience  has  also  shown  that  when  the  signal  is  scintillating  violently  the 
tracking  filter  will  not  stay  locked.  For  these  reasons,  and  the  difficult  align- 
ment procedure,  it  appears  that  frequency  tracking  is  not  the  solution  for  low-cost 
satellite  monitoring  receivers,  although,  in  most  cases,  it  is  the  only  solution  at 
present. 


2.  PHASE  MEASURING  EQUIPMENT 

Fortunately  there  are  satellites  transmitting  high  and  low  UHF  signals  that 
are  coherent;  that  is,  both  are  multiples  of  the  same  crystal  frequency  and,  as  a 
consequence,  lend  themselves  favorably  to  a simpler,  much  more  dependable 
and  less  expensive  system  for  their  detection  and  processing.  This  is  a cross - 
correlation  technique  designed  at  the  AFCRL  Ionospheric  Research  Laboratory’, 
Bedford,  MA  to  measure  differential  phase  between  a 40 -MHz  signal  and  a 
360-MHz  signal  transmitted  from  the  ATS-6  satellite.  If  the  use  of  electronic 
correlation  is  not  restricted  because  of  the  required  response  time,  it  is  a power- 
ful tool  capable  of  extracting  minute  but  meaningful  information  from  an  excep- 
tionally noisy  environment.  The  bandwidth  determined  by  the  integration  time 
constant  can  be  made  extremely  narrow.  In  this  system  it  is  0.  ,5  Hz  at  -.1  dB, 
this  being  adequate  for  the  fastest  360°  phase  shift  encountered. 

Inspection  of  the  cross-correlation  formula,  i (r>  - l/T  j f{T)  g(t  - r)  dt  shows 
that  when  two  time-varying  voltages  f(t)  and  g(t  - t)  are  the  same  frequency  and 
the  delay  time  T is  zero,  the  signals  will  be  in  phase  and  the  output  voltage,  ® (t) 
will  be  a maximum  positive  dc  voltage.  As  T varies,  the  phase  between  the 
voltages  increases.  At  90°  the  average  correlated  output  is  zero.  When  T is  such 
that  the  phase  difference  is  180°,  the  correlated  output  is  a maximum  negative  dc 
voltage.  If  T is  varied  linearly  in  incremental  steps,  using  a tape  loop,  for 
instance,  to  provide  g(t  -r),  the  correlated  integrated  output  is  sinusoidal  with 
each  incremental  step's  worth  of  integrated  data  being  produced  for  each  cycle  of 
the  tape  loop. 

The  use  of  cross-correlation  to  detect  and  process  ionospheric  phase  data 
provides  a compatible  marriage  between  the  processing  technique  and  the 
phenomena  producing  the  ionospheric  anomalies  for  the  following  reasons: 

(1)  The  transmitted  frequencies  are  made  coherent  in  the  satellite,  therefore, 
frequency  locking  equipment  is  not  required  for  this  equipment. 

(2)  The  variable  delay  (t)  does  not  require  any  additional  equipment  as  it  is 
generated  naturally  in  the  ionosphere  by  the  diurnal  variation  of  the  total  electron 
content,  with  its  inherent  refractive  effects  changing  the  phase  velocity  of  the 
40-MHz  signal.  T is  also  sensitive  to  phase  scintillations  and  faster  changes  in 
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F'igure  1.  40-  and  360-MHz  Differential  Phase  System  Block  Diagram 

the  total  electron  content  (TEC),  such  as  those  caused  by  magnetic  storms  and 
ionospheric  winds;  also  traveling  ionospheric  disturbances,  that  is,  acoustic- 
gravity  waves  that  modulate  the  TEC  profile  as  they  propagate  through  the 
ionosphere.  ^ 

(3)  Because  the  system  rejects  random  noise  and  uncorrelated  signals,  the 
products  of  which  average  out  to  zero,  the  normal  receiver  requirements  are 
minimal,  and  it  is  not  necessary  to  employ  a complex  phase  locked  communications 
receiver. 

(4)  The  circuitry  that  provides  a 0.  5 Hz  bandwidth  also  produces  a quasi- 
sinusoidal  integrated  output  voltage  when  the  TEC  is  increasing  or  decreasing  at 
a linear  rate.  By  coincidence,  and  very  conveniently,  this  satisfies  the  first 
requirement  for  generating  a voltage  proportional  to  phase  for  making  a strip 
chart  recording. 

A block  diagram  (Figure  1)  has  been  provided  to  assist  in  the  explanation  of 
the  complete  system.  The  receiver  portion  consists  of  blocks  1 through  7.  Blocks 
1 and  4 are  standard  preamps  with  a 2 dB  noise  figure  and  a 28  dB  gain.  The 
local  oscillator,  block  7,  is  tuned  2 KHz  higher  than  the  40-MHz  satellite  signal; 
thus,  when  mixed  with  the  satellite  signal,  2 KHz  becomes  the  difference  frequency 
for  further  processing. 

1.  Yeh,  K.C,  and  Liu,  C.  H.  (1974)  Acoustic -Gravity  waves  in  the  upper 
atmosphere.  Rev.  Geophys.  Space  Phys.  12(No.  2):193. 
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Figure  2.  Square  Law  Spacing  of  I'haso 
Cycles  Indicating  a I, inear  Separation 
Kate  of  the  3fiO-MHz  and  Kefracted 
40 -MHz  Rays 


Uiien  the  local  oscillator  is  mixed  with  the  3G0-MHz  satellite  signal  it  is 
actually  the  Oth  harmonic  that  is  used,  thus  making  the  2 KHz  difference  frequency 
Ifl  KHz  for  further  processing  at  the  mixer  output. 

At  this  point  the  system  has  produced  a 2-KHz  signal  carry  ing  the  40-MHz 
phase  information  and  an  18-KHz  signal  to  be  used  as  the  standard,  since  refrac- 
tive effects  on  the  360-MHz  frequency  responsible  for  it  are  minimal. 

The  two  frequencies  have  remained  coherent  through  the  heterodyning  process, 
because  the  same  L.  O.  has  been  employed  for  both;  however,  for  cross- 
correlation  purposes  they  should  both  be  the  same  frequency.  To  accomplish  this 
the  18  KHz  could  be  divided  down  to  2 KHz  or  the  2 KHz  could  be  multiplied  up  to 
18  KHz.  The  latter  method  was  chosen  for  two  reasons.  There  is  more  stability 
in  frequency  multiplication  than  there  is  in  frequency  division,  and  it  is  possible 
to  make  a chart  record  that  is  nine  times  more  sensitive,  nine  being  the  multiplica- 
tion factor  to  bring  2 KHz  up  to  18  KHz.  Therefore,  each  ®(t)  cycle  out  of  the 
correlator  represents  40’  instead  of  360°  of  the  40-MHz  phase  change.  See 
Figure  2. 

The  increased  sensitivity  is  valuable  when  measuring  the  time  of  occurrence 
of  TUI's  (Figure  2)  as  the  phase  reversals  can  be  measured  accurately  to  within 
a few  degrees.  Its  potential  value  is  the  accuracy  to  be  achieved  in  measuring 
TUI  direction  and  speed,  if  three  systems  are  put  into  operation  in  triangulation. 

For  strip  chart  recording  of  phase  and  phase  scintillation,  a sinusoidal 
function  is  not  adequate  as  it  is  subject  to  ambiguity,  that  is,  one  cannot 
tell  if  the  phase  is  increasing  or  decreasing.  If  the  system  includes  two  cor- 
relators, one  generating  a sinusoidal  function  and  the  other  generating  a cos 
function,  as  this  system  does  (blocks  10  through  13),  the  ambiguity  can  be 
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eliniinatod  thfougli  visual  interpretation  of  the  record.  However,  this  is  a tedious 
operation,  so  it  is  better  to  extend  the  process  further  and  include  a coordinate 
converter  (l)locks  12,  14,  15,  I G,  17)  that  puts  the  data  back  in  the  form  of  two 
coherent  ac  signals  still  retaining  the  proper  phase  relationship  that  can  operate 
a phase  ramp  generator.  The  frequency  can  be  whatever  one  chooses.  This 
system  uses  GO  Hz  because  of  future  intentions  to  drive  a clock  motor.  It  may 
appear  at  this  point  that  nothing  has  been  gained;  we  started  with  two  18-KHz 
signals,  one  changing  phase  relative  to  the  other,  and  we  end  up  with  two  GO-Hz 
signals,  one  changing  phase  relative  to  the  other.  The  big  difference,  however, 
is  that  the  unprocessed  18-KHz  signals  have  so  much  phase  jitter  and  extraneous 
noise  that  the  data  from  a phase  detector  operating  on  these  signals  would  be 
worthless.  On  the  other  hand,  the  two  GO-Hz  signals  are  relatively  noise-free 
having  gone  through  the  correlation  process  with  the  0.  5 Hz  bandwidth. 

Mathematical  expressions  describing  the  more  complex  portions  of  the  system 
are  as  follows:  Placing  a 90^  phase  shift  on  one  of  the  18-KHz  signals  and  simul- 
taneously running  it  through  a second  correlator  produces  a second  correlated 
signal,  thus  providing  two  output  voltages,  sin  0(t)  = 1/T  f(t)  g(t  -t)  dt  and 
cos  j(r)  = 1/T  *^f(t)  g(t  -T)dt. 

Jo 

When  the  phase  between  the  two  input  signals  is  shifting  at  a constant  rate, 
these  two  functions  can  be  considered  to  be  very  low  sub  audio  frequencies,  one 
90°  out  of  phase  with  the  other.  One  cycle  per  second  is  about  the  highest  fre- 
quency encountered  so  far.  Therefore,  treating  them  as  such  yields  the  following 
expressions  for  the  output  of  the  modulators:  box  14,  (sin  xjt)  (sin  i.2t)  and  box 
16,  (cos  1 it)  (cos  X 2^)-  Summing  the  two  outputs  (box  17)  yields  the  trigono- 
metric relation,  cos  (i'  it  - i 2*^-  Amplifying  this  signal  into  saturation  obliterates 
the  cos  function  and  the  output  is  then  tCi:^  - x ^).  Thus  the  output  of  the  summer 
and  amplifier  when  'r;^  = K = 60  Hz  and  j:  ^ ~ 1 Hz  is  a 59-Hz  square  wave  which  is 
just  another  way  of  saying  a 60-Hz  constant  amplitude  voltage  is  changing  phase 
at  the  rate  of  one  cycle  per  second  relative  to  the  60-Hz  standard. 

The  multiplication  by  l/T  is  performed  to  normalize  the  data,  that  is, 
to  eliminate  signal  time  duration  as  a factor  in  determining  the  amplitude  of  the 
integrated  output. 

Operating  a conventional  phase  detector  off  of  the  two  GO-Hz  signals  produces 
a phase  ramp  for  each  40°  phase  shift  between  the  two  rf  signals  at  the  antenna, 
that  is,  360°/ 9 as  explained  earlier  in  the  text. 

It  was  decided,  however,  not  to  phase  detect  here,  but  to  extend  the  range  to 
1280°  per  ramp  to  help  distinguish  phase  scintillations  from  phase  ramps  and  also 
to  reduce  the  chart  speed  for  more  economical  use  of  chart  paper.  To  accomplish 
this  both  60-Hz  signals  are  fed  into  32  to  1 digital  counters,  (boxes  19,  20).  The 
counter  output,  operating  from  the  60-Hz  oscillator,  supplies  a pulse  to  repeatedly 


9 


W. 


I'lguri’  Itiflucnco  of  Triivclinf;  lono- 
sphi-ric  l)isturl)atK-i'  on  l^lectron  (!(.)ntcnt 
,in<l  HI'  I’roijagalion  1 sing  (ieomotrx  at 
li(tuator  for  Simplification 


start  a very  linear  ramp  generator  (block  21).  This  is  a sawtooth  voltage  with  : 
repetition  rate  of  approximately  0.  f)  cycles  per  second,  that  is,  .12/00  Hx. 

The  ramp  generator  is  fed  to  the  phase  detector  where  it  is  sampled  h\  the 
instantaneous  pulse  from  the  22  to  1 signal  counter.  The  value  sampled  is  placeil 
in  storage  on  a capacitor  and  held  constant  until  the  next  pulse  2 seconds  later. 

If  the  phase  is  increasing,  the  sampled  values  follow  the  ramp  up  in  the  positive 
direction.  If  the  jihase  is  decreasing,  the  reverse  is  true.  I'or  example,  if  it 
takes  12  seconds  for  the  signal  Co  coniplete  each  cycle  of  the  2-second  ramp  in 
incremental  stored  steps,  the  output  of  the  phase  detector  is  a il2-second  ramp 
that  reverses  direction  when  the  signal  reverses  phase.  See  f igure  2. 

Karlier  in  the  text  it  was  mentioned  that  there  was  a reason  for  using  GO  I lx 
in  the  coordinate  converter.  The  proposal  is  to  develop  a complete  diurnal  curve 
of  the  total  electron  content  using  electromechanical  techni(^ues.  If  the  two  GO-llx. 
voltages  (blocks  12  and  17)  are  used  to  drive  clock  motors,  the  output  shaft  of  the 
signal  GU  llz  will  rotate  slightly  faster  or  slower  than  the  standard  fiO-lIz  motor, 
depending  on  how  fast  the  phase  is  increasing  or  decreasing.  Hy  using  an  appro- 
priate gear  reduction  in  combination  with  a mechanical  differential,  coupled  to  a 
linear  potentiometer,  it  will  he  possible  to  develop  a variable  ric  voltage  repre- 
senting a diurnal  plot  of  the  electron  content,  tliere  being  a direct  correlation 
between  phase  shift  and  changing;  Tlit:.  The  gear  reduction  will  have  to  be  suf- 
ficient to  keep  the  rotation  of  the  potentiometer  arm  within  the  mechanical  limits 
of  the  potentiometer.  This  will  ensure  a smooth  curve  without  interruption 
caused  by  potentiometer  discontinuity. 
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! 'it;urc  t.  St;itisl it  at  CnuiiKirisoii  of  Ami->litiidc’  and 
I’hasc  Scintillation 


1.  ( <)\U’\KIN(,  Wtl’l.m  l)K  SCIMIM.MION  WITtl  IMIASl.  St.lM  IM.M  ION 

I'lio  equivalent  tliin  phase  screen"  is  a popular  concept  used  by  iono- 
spheric scientists  for  explaining,  in  niatheniatical  terms,  amplitude  scintillations 
on  rf  signals  propagated  through  the  ionosphere.  This  is  a screen  of  negligible 
thickness  that  would  produce  the  same  phase  variations  on  the  rf  signal  as  the 
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actual  irregular  region.  “ Tlie  assumption  is  tliat  the  signal  leaves  the  thin 
screen,  fluctuating  in  phase  only  hut,  as  the  wave  propagates  beyond  the  screen, 
amplitude  fluctuations  develop  from  the  summation  of  the  deviated  phase  signals. 

Kecent  measurements  made  on  ATS-fi  geostationary  satellite  signals  com- 
paring 10-MMz  phase  with  amplitude  sliow  that  application  of  the  equivalent  thin 
screen  theorv  is  valid  for  only  a portion  of  the  data.  During  a 13-day  active 
[tenod  of  ionospheric  scintillation  in  \uvember  1!174,  data  was  computerir.ed  to 
make  the  above  comparisons.  I'iie  results  show  that  during  the  periods  of 
amplitude  scintillation,  phase  scintillation  existed  only  33  percent  of  the  time 
(I'igtire  1).  I'hus,  if  the  amplitude  scintillation  at  the  antenna  was  not  the  result 
of  ionosphi'ric  phase  scintillation  during  47  percent  of  the  time,  another  theory  is 
required  to  explain  the  results.  To  conform  to  the  data  as  presented  on  the  strip 
charts,  the  responsiiile  mechanism  must  allow  both  enhancement  and  attenuation 
of  the  rf  signal.  I’ast  investigators  studying  large-scale  irregularities  have 

2.  I'.rigg.s,  H.  n.  (1  ;ir>r.)  Ilrief  review  of  scintillation  studies,  Hadio  Sci. 

l(\o.  1 01:1103. 


1 I 


1 

40-WMi 

PHASE 


I'igure  5.  AmplitutU*  ami  Phase  Scintillation  Occurring 
T ogether 


suggested  marked  focusing  effects  as  the  cause  of  scintillation.  This  could  also 
explain  the  remaining  47  percent  reported  here.  Concurrent  cycles  of  the  rf 
propagation  could  be  focused  on  the  antenna  through  lens-like  structures  of  the 
irregularities,  eliminating  the  phase  deviation  requirement  dictated  by  the  dif- 
fraction theory.  Also,  when  a moving  ionospheric  structure  focuses  an  increase 
in  energy  on  the  antenna,  its  refractive  area  that  supplied  the  enhancement  must 
also  reduce  the  energy  when  it  is  appropriately  positioned,  resulting  in  moving 
irregularities  that  are  equivalent  to  alternating  convex  and  concave  lenses. 

In  a similar  manner,  moving  electron  gradients  normal  to  tlie  propagation 
path  could  deflect  concurrent  rays  that  are  close  to  grazing  incidence  resulting  in 
alternating  enhancement  and  fading  of  the  signal  at  the  antenna.  ^ An  exatnple  of 
amplitude  scintillation  accompanied  by  phase  scintillation  is  shown  in  Figure  7; 
amplitude  scintillation  without  the  presence  of  phase  scintillation  is  shown  in 
Figure  6. 

,A  graph  (Figure  7)  was  constructed  showing  the  percentage  of  time  the  phase 
angle  ; was  in  excess  of  x degrees  for  the  13-clay  period.  The  sniootlmess  of  the 
curve  suggests  that  sufficient  data  was  analyzed  to  describe  conditicins  for  a higl: 
percentage  of  ionospheric  scintillation. 

Although  the  shape  of  the  curve  appears  to  indicate  a near  Gaussian  distribu- 
tion of  irregularities  with  respect  to  the  line  of  sight  between  the  satellite  and  the 
receiver,  logically  there  is  no  reason  to  believe  there  would  he  such  a distribution 
at  the  specific  region  where  the  rf  ray  penetrates  the  ionosphere.  It  is  more 
likely  that,  statistically,  there  is  a linear  distribution  being  affected  li\  a dimin- 
ishing power  law. 

3.  .Slack,  F.  F.  (1!)71)  Quasiperiodic  scintillation  in  the  ionosphere,  J.  Atmos. 

Terr.  I’hys.  34:927. 
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Figure  6.  Amplitude  Scintillation  Without  Phase 
Scintillation 


Figure  7.  Phase  Scintillation  in 
Degrees 


The  Fresnel  formula,  vt  = ^'2nX r states  that  when  the  lateral  distance  in  the 
ionosphere  lietween  the  line  of  sight  ray  and  an  interfering  ray,  or  the  equivalent 
radius  resulting  from  diffracted  or  scattered  rays  has  changed  vt,  it  produces  a 
difference  in  the  arrival  time  required  for  ni  at  the  receiver  distance  r.  If 
j^/2'*  is  substituted  for  n,  the  equation  can  be  solved  for  phase  change  instead 
of  integers  of  wavelengths.  Thus  with  rearrangement  and  substitution,  the  formula 
liecomes 


I'rnm  this  it  is  :<ppar(>nt  that  the  change  in  pliase  resulting  from  a linear  distribu- 
tion IS  proportional  to  the  square  of  the  ray  sepai-ation  distance  vd.  Looking  again 
at  tile  grapli  in  tlie  region  of  > I rad,  the  curve  closely  follows  a square  law 
indicating  I'resnel  action  on  a linear  distribution  of  irregularity  sizes.  If  such  is 
the  case,  it  supports  tlie  contention  made  by  several  investigators  that  some  iono- 
st'iieric  scintillation  is  produced  liy  either  a superposition  or  a distortion  of  quasi- 
periodic  ('\eiits.  Tliese  arc  the  irregularities  which,  when  isolated,  develop  a 
characteristic  ringing"  I'resnel  pattern  as  displayed  on  ttie  amplitude  channel  of 
a chart  recorder.  *'  ^ They  fall  in  the  region  » 1 rad,  and  the  curve  would 
liave  to  lie  extended  significantly  to  include  the  uncontaminated,  isolated  examples. 
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■file  region  of  T < 1 rad  is  where,  according  to  Ilewish,  and  generally 
accepted  li\  otliers,  the  projected  irregularities  have  the  same  scale  size  when 
mea.stired  on  the  ground  as  would  be  the  case  if  measured  at  the  ionospheric  lev'el. 
Mere  the  curve  is  quite  linear,  as  one  might  expect  with  a linear  distribution  of 
increasing  scale  sizes,  because  the  cliange  in  refractive  index,  and  thus  the  cliange 
in  jihase,  is  directly  proportional  to  th<>  electron  content  in  the  irregularity. 

.Aarons  el  al'  in  referring  to  the  two  regiiiK-s  as  weak  scattering  wliere  t •'  1 rad 
and  strong  scattering  wliere  1 rad  shows  that  there  are  periods  when  the  two 

regimes  exist  simultaneously.  The  gra|>h  supports  their  findings  by  displaying  a 
gradual  transition  from  linear  to  a eluiracteri Stic  l-'resnel  curve  in  the  region 
between  'i7  and  ''0°. 

■['he  reeordefl  data  does  not  diff(*rentiato  between  positive  and  negative  phase 
deviations,  so  they  till  have  been  lumped  on  one  side  in  the  graph.  Since  7^  is  the 
result  of  a ehanging  refractive  index  about  a mean  level,  a complete  graph  sliould 
include  a mirror  image  of  itself  in  the  m'galive  7 diri'ction  to  more  accurately 
(lortray  the  process  that  produces  the  I'hase  scintillations.  This  is  taken  itito 
consideration  in  explaining  the  non-linearity  of  the  curve  close  to  ' 0.  It  is 
similar  to  ttie  zero  beat  that  apiiears  on  a traveling  ionosplierie  disturbance 

4.  Likins,  T..!.  and  Slack,  F.  L.  (1 9G9)  ( ibservations  of  traveling  ionospheric 
disturliancos  using  stationary  satellites,  .1.  Atmos.  Terr.  Phys.  :ll;421. 

a.  rithcridge,  .1.  F.  (11)71)  Tlie  diffraction  of  satellite  signals  by  isolated  iono- 
spheric irregularities,  .1.  Atmos.  Terr.  Phys.  4.1:47. 

fi.  Kelleher,  11.  F.  and  Martin,  P.  (197.'))  Fresnel-type  fading  on  satellite 
records  at  low  latitudes,  .1,  Atmos.  Terr.  Phys.  47:1109, 

7.  Hewish,  A.  (lO.'il)  The  diffraction  of  radio  waves  in  passing  through  a phase- 

changing ionosphere,  Proc,  H,  Soc.  ■\20'):81. 

8.  Aarons,  .J.,  Castelli,  J.  P. , and  Kidd,  VV.  (19G2)  Cygnus  rise  and  set 

measurements  at  VHP  and  UHF  radio  frequencies,  Proc.  Internat,  Conf. 
lonos.,  252. 
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record  (I'lKure  3)  when  the  phase  reverses  direction,  and  tlie  angular  refraction 
goes  to  zero. 

In  comparing  the  information  contained  in  Figures  4 and  7,  and  using  a spate 
of  speculation  with  at  least  a modicum  of  logic,  it  appears  that  there  is  a relation 
betweeti  : scintillation  in  the  hours  between  sunset  and  sunrise  and  < 1 rad. 

•Since  ttie  electron  content  is  decaying  during  this  period,  and  the  electrons  would 
be  clinging  to  tlieir  magnetic  field  aligned  positions,  electron  gradients  would 
develop  normal  to  tlie  field.  With  the  ionosphere  in  motion,  conditions  now  e.icist 
for  the  occurrence  of  c scintillations,  initially  falling  in  the  regime  < 1 rad, 
because  adjacent  areas  would  produce  only  slight  differences  in  refraction.  As 
the  decay  cycle  continues,  electron  deficient  holes  would  likely  develop  because 
of  the  remaining  electrons'  propensity  for  field  alignment,  causing  the  medium 
to  become  more  irregular.  At  this  time  there  would  be  a tendency  towards  focusing 
meclianisms  or  grazing  reflection  surfaces  more  suitaljle  for  t scintillations  in 
tile  regime  > 1 rad.  Hriggs  in  writing  about  "the  equivalent  thin  phase  screen" 
and  its  application  to  scintillation  amplitude,  summarized  by  saydng,  "This  method 
will  give  useful  results  if  < 1 rad.  When  is  large  there  is  evidence  of  focus- 
ing. " .\ote  in  Figure  4 the  negative  correlation  between  j scintillation  and  ampli- 
tude scintillation  during  the  daylight  hours.  Then  compare  this  with  Figure  7 
showing  that  the  rate  of  occurrence  of  j scintillation  is  practically  zero  for 
240^  . Thus  the  data  indicates  that  during  the  daylight  hours  the  increased 
electron  density  requires  that  the  irregular  structure  instead  of  behaving  like  holes 
must  now  behave  as  cylinders  with  essentially  all  amplitude  scintillations  being 
caused  by  focusing  or  reflection  surfaces,  that  is,  i scintillation  is  almost 
nonexistent,  and  when  it  does  occur  it  is  likely  to  be  in  the  regime  ® >1  rad. 

It  must  be  emphasized  that  the  above  explanation  has  been  based  on  13  active 
days  of  scintillation  in  November  and  its  application  to  scintillation  in  general  is 
a matter  of  conjecture. 

4.  THK  TR.WF.I.INC  lOlNO.sPHKRIC  Di.STl'RBANCK 

One  of  the  most  outstanding  features  displayed  by  the  phase  data  is  the 
influence  that  traveling  ionospheric  disturbances  (TID)  have  on  the  phase  of  the 
40 -MHz  signal.  A TID  is  generally  thought  to  be  an  acoustic -gravity  wave 
propagating  through  the  vast  expanse  of  the  atmosphere  including  the  ionosphere. 

Interpretation  of  the  data  (Figure  2)  indicates  that  the  medium  is  being  dif- 
fused and  compressed  in  a cyclic  manner  that  changes  its  refractive  index  in  the 
same  fashion.  I'urther,  the  data  shows  that  the  alternating  diffusion  and  compres- 
sion is  in  a direction  that  crosses  the  satellite  receiver  propagation  path.  To 
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Figure  8.  TKC  Slant  Column 


1 1 >.c  li  tins  conclusion  consider  tlic  diurnal  curve  of  the  total  electron  content 

(i  u uie  iO.  This  IS  the  sum  of  all  the  phase  ramps  into  one  continuous  24-hour 

I I I \ ( , Since  there  is  a direct  correlation  between  TEC  and  phase,  the  ordinate 

2 

I I tin  curve  ciin  lie  calibrated  in  tiie  nunilior  of  electrons  per  m column  between 
n ( s;  lellite  and  tlie  receiver.  I ’oliirimeter  data  has  lieen  used  this  way  for  a 
I iin  1 of  \cars  but  uses  the  reltitionslhp  ijetween  TEC  and  Faraday  rotation  and 
1-  oiiK  uliout  10  )ieroent  as  sensitive  to  the  changes  in  the  medium.  Also  the 
!•  .1.0  iiiietcr  measures  tlie  I arada>  rotation  only  to  the  range  still  influenced  liy 
b ( urlb's  magnetic  field.  Hecause  of  this  difference  the  electron  content,  as 
I < ,i.-Mir<  d b,\  the  polarlmeter,  should  never  exceed  the  differential  phase  TEC 
1.0  . surements.  Since  the  figure  contradicts  this  fact,  it  appears  that  the  wrong 
I . ij  e value  was  selected  in  the  polarimoter  analysis. 

1 lie  figure  shows  that  the  Til)  produce  equal  deviations  about  the  median  of 
1 ..  ( urve.  riiis  indicates  tlial  (lic'-e  has  been  no  net  gain  or  loss  in  the  total 
1 I.  I I /on  l ontent,  o/d_\  tliat  the  electrons  have  been  diffused  away  from  t ie  sat- 
,,  , . ui  r I'lopagation  path  ami  then  compressed  back  again.  If  the  compression 
;.'.ii  diftusion  were  in  the  direction  of  the  rf  path,  the  data  would  show  no  change 
,1.  t'.(  total  electron  content  column  when  Til)  are  present. 

.Since  the  ion  gv rofrequenev  i.s  much  higher  than  both  the  ion  neutral  collision 
fi  I iiueni  \ and  the  wave  frequency,  the  charged  particles  can  move  only  along  the 
inaciietic  field  lines.  * Thus  ionospheric  oscillations  are  induced  more  easily  by 
Ik  Id  aligned  perturbations.  Combining  this  information  with  the  above  presents 
the  hypothesis  that  the  extent  of  Til)  influence  on  the  rf  propagation  path  is 
latiturle  dependent,  being  maximum  at  the  equator  where  the  satellite-ground 
propagation  path  is  normal  to  the  magnetic  field  lines,  and  minimum  in  the 
hoi  them  latitudes  where  the  ray  path  and  field  lines  are  more  closely  parallel. 
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1 , i f cdut  bc.  I ncludes  otlu  r factors  that  may  be  more  overpowering,  such  as 

^ TID  having  spent  its  energy  by  the  time  it  reaches  the  equator. 

Ttio  ATS-G  satellite  elevation  angle  at  Sagamore  Hill,  Hamilton,  Massachusetts 
IS  GG'’.  The  magnetic  dip  angle  is  73°.  Therefore  the  rf  propagation  path  crosses 
llic  magnetic  field  lines  at  an  angle  of  37°.  It  follows  that  the  magnitude  of  the 
electron  flow,  normal  to  the  propagation  path,  is  sin  37°  k the  magnitude  of  the 
electron  motion  moving  back  and  forth  along  the  field  lines  as  they  respond  to  the 
TID  wave  motion.  This  compares  with  sin  90°  at  the  equator. 

It  is  also  important  to  consider  the  TID  direction  component  normal  to  the  rf 
propagation  patli,  so  that  the  refractive  process  taking  place  can  be  understood, 
since  this  is  what  the  receiver  and  phase  detector  system  is  measuring.  If  the 
crest  of  the  TID  wave  represents  the  maximum  electron  density  and  the  trough 
till'  minimum  (Figure  2),  then  one  can  visualize  this  as  introducing  ionospheric 
gi  adients  in  the  direction  of  the  TID  motion  with  maxima  and  minima  representing 
positions  wliere  the  rate  of  change  of  electron  content  goes  to  zero.  Zero  angular 
M.  fraction  (no  beniling  of  the  rf  ray)  would  appear  when  these  points  are  centered 
i"i  tiiv  ray;  more  specifically,  when  components  normal  to  the  ray  have  reached 
t'  I imdition,  that  is,  the  integrated  column  of  electrons  are  equal  on  opposite 
b'  h s of  tlu-  ray  m the  incident  plane.  It  follows  then  that  maximum  refraction 
I I 1 111  ,s  when  the  rate  of  cliange  of  electrons  is  greatest  during  the  steepest  part 

0 t)  ( ’I'll)  uaio.  'Diis  is  Just  what  the  recorded  data  shows,  except  that  these 

1 I. mill  ions  tiro  sliglitly  offset  l>v  the  ever  present  diurnal  change  in  total  electron 
content. 

I he  earlit'f  section  on  the  equipment  shows  how  these  physical  interactions 
t II  1 1 analatod  into  electrical  measurements  that  are  displayed  on  a chart  recorder. 
^ In  tier  understanding  is  provided  by  studying  the  diagram  (Figure  3).  The  dia- 
t leiii  represettts  conditions  in  the  ionosphere  where  the  satellite  is  directly  over- 
hind;  the  magnetic  field  lines  are  at  right  angles  to  the  rf  propagation  path,  and 
till  TH)  direction  is  also  normal  to  the  rf  propagation  path.  The  figure  shows 
M h.'d  iiappens  to  the  rf  ray  as  a TID  influenced  electron  gradient  moves  across  it. 

.\'p  bending  occurs  at  tlie  maximum  and  minimum  positions  where  there  is  zero 
gradient.  Maximum  liending  occurs  during  the  steepest  portion  of  the  gradient. 
Since  the  refractive  index  is  less  than  one,  the  bending  is  always  towards  the 
greater  density. 

As  the  gradient  passes  across  the  ray,  the  phase  velocity  completes  a cycle 
of  increasing,  decreasing,  and  reversing  direction,  then  increasing  and  decreas- 
ing in  the  reverse  direction.  One  side  is  a mirror  image  of  the  other  except  for 
tlio  imlialance  caused  t)v  the  diurnal  effect,  or  some  other  natural  perturbation. 

This  action  repeats  for  each  passing  TID  cycle  that  intercepts  the  rf  ray  and  is 
illustrated  on  the  strip  chart  (Figure  2). 
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There  have  been  pen  recordings  showing  examples  of  TID  electron  displace- 
ments up  to  20  percent  of  the  ambient  total  electron  content  at  the  time  of  the 
disturbance.  Close  inspection  of  the  phase  data  40°  cycles  and  1280°  phase 
ramps  shows  that  the  refraction  generally  is  causing  the  40-  and  3G0-MHz  rays 
to  separate  and  close  at  a constant  velocity,  which,  in  turn,  causes  typical  doppler 
spacing  of  the  40°  cycles.  The  cycles  follow  the  square  law  spacing  for  a con- 
stant velocity  fixed  frequency  moving  at  right  angles  to  the  propagation  path.  A 
mathematical  solution  shows  that  doubling  the  lateral  distance  the  ra\  travels 
produces  four  times  the  number  of  cycles.  Note,  on  Figure  2,  that  doubling  the 
2-C3'cle  distance  produces  8 cycles,  and  doubling  the  4-cycle  distance  produces  IG 
cycles.  It  is  assumed  that  the  dominant  action  is  on  the  40-MHz  ray,  the  coherent 
3G0-MHZ  being  high  enough  to  be  almost  free  of  the  refractive  effects  and  is  thus 
used  as  the  standard  reference  for  the  phase  measurement. 

It  should  be  emphasized  that  the  number  of  cycles  between  crest  and  trough 
per  unit  time  does  not  measure  the  velocity  of  the  TID,  but  relates  to  the  rate  of 
separation  and  closing  of  the  40-  and  3G0-MHz  rays  normal  to  the  propagation  path. 
The  number  of  cycles  are  also  a measure  of  the  TID  amplitude;  that  is,  the 
bigger  the  TID,  the  greater  the  influence  on  the  electron  density  wliich,  in  turn, 
produces  a larger  change  in  refractive  index  and,  thus,  an  increase  in  phase 
velocity,  as  illustrated  on  the  chart,  by  recording  a greater  number  of  cycles  per 
unit  time. 

Doppler  theory  adequately  explains  the  form  of  the  data  as  displayed  on  the 
chart  record;  however,  many  of  the  parameters  that  influence  space  communica- 
tions can  be  gleaned  from  the  data  through  the  application  of  ray  optics  analysis, 
which  then  permits  the  use  of  Snell's  law  and  the  Appleton  Hartree  formula.  A 
list  of  these  parameters  in  addition  to  totad  electron  content  (TKC)  are: 

(1)  Changes  in  electron  density,  assuming  height 

(2)  Slab  thickness,  assuming  heignt 

(3)  A continuous  measurement  of  phase  velocity 

(4)  A continuous  measurement  of  rate  of  change  of  TKC 

(5)  A continuous  measurement  ot  refractive  index 

(6)  A continuous  measurement  of  angle  of  arrival 

A simple  analysis  can  be  performed  on  the  recorded  data,  on  those  sections 
where  the  chart  record  clearly  indicates  that  the  phase  is  shifting  in  accordance 
with  a linear  motion  normal  to  the  rays;  that  is,  the  40°  cycles  follow  the 
square  law  spacing  mentioned  earlier.  There  are,  generally,  several  such 
sections  to  be  found  each  day. 

Ray  optics  analysis  is  valid  if  the  following  conditions  are  met;  (1)  the 
incident  angle  is  sufficiently  large;  (2)  the  medium  is  lossless:  (3)  the  phase  is 
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I’igure  9.  The  Geometry  for  the  Case 
Where  the  Satellite  Is  Oirectly  (X'erhead 
and  the  TID  Direction  Is  Normal  to  the 
rf  Hay 


measured  at  a distance  sufficiently  far  from  the  inhomogeneous  medium;  and 
(4)  the  media  vary  slowly  compared  to  a wavelength.'  These  conditions  exist 
for  the  data  presented.  The  first  condition  exists  because  gradients  normal  to 
the  ray  patli  caused  by  the  TID  arc  considered  to  be  the  dominant  influence  on  the 
rf  ray:  tiierefore,  tiie  incident  angle  ?.  is  measured  from  a line  normal  to  the 
gradient.  Figure  9 shows  that  7.  is  always  very  large  going  to  90^  at  both  tlie 
crest  and  trough  of  tlie  TlfJ.  Since  total  reHection  occurs  in  all  dielectrics  at 
the  grazing  incident  angle  (90’ ),  the  principle  of  reversitiility ' ^ can  lie  used  to 
explain  tlie  phase  reversal  at  botli  the  crest  and  trough.  The  principle  of  rever- 
sibility states  that  whenever  a ray  is  totally  reflected  a phase  change  of  180' 
occurs;  also,  if  the  phase  reversal  docs  not  take  place  when  the  wave  passes  from 
the  lower  to  the  higher  refractive  index,  it  will  occur  when  the  wave  passes  from 
the  higher  to  the  lower  refractive  index.  This  explains,  through  ray  optics,  the 
phase  reversals  appearing  on  the  chart  rccoivls  at  both  the  ITD  crest  and  trougli 
riius,  ray  optics  tfieory  provides  the  same  answer  given  by  the  y)hase  doppler 
nroecss  which  shows  the  frequency  or  phase  to  be  increasing  when  range  of  the  rf 
source  is  decreasing,  and  decreasing  when  the  range  is  increasing,  with  a phase 
reversal  at  or  near  the  point  of  closest  approach.  Figure  9 depicts  the  geometry 
for  the  case  where  the  satellite  is  directly  overhead  and  the  TID  direction  is 
normal  to  the  rf  ray. 


9.  Papa,  R..T. , private  communication. 

10,  Jenkins,  F,A.  and  White,  H.  K.  (1937)  Fundamentals  of  Optics,  McGraw- 
Hill,  New  York. 
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Solving  the  geometr>-  yields  some  important  relationships: 


iyt)^ 

2ht 


and 


m = s 2n> r.  (2) 

This  is  recognized  as  the  Fresnel  formula  when  n - 1,  and  states  that  the  ray 
separation  or  closing  must  have  traveled  t/t  in  order  to  change  the  phase  360’  at 
the  distance  r. 

If  the  incident  angle  is  measured  from  the  path  ui,  a line  normal  to  the 
TID  gradient,  the  relationship  in/r  = cos  j.  is  important,  because  it  accom- 
modates the  introduction  of  an  equation  derived  from  Snell's  law,  cos  - f^^/ f 

where  f.,  = plasma  frequency  and  f = the  propagation  frequency.  From  the 

2 2 

Appleton -Hartree  equation:  u^  = l-  x = l-  (f^,/f)^  = 1 - KN7.2,  where  u = the 

2 2 ^ 3 

refractive  index,  K = (e  !\r  c^m)  = 80.6,  f is  in  Hz,  N = number  of  electrons/ m 

at  highest  density  point.  Combining  these  equations  gives: 
f^.  = ./80.  6N 
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2nXf^ 

■ 80.  ^ 


or  N 


. nXf^ 

■ 40.  3r  • 


(3) 


The  last  expression  states  that  if  the  range  to  the  point  of  highest  electron 
density  is  known,  the  change  in  density  caused  by  the  TID  cam  be  calculated. 
However,  it  must  be  remembered  that  the  equation  represents  conditions  wtiet: 
the  40-  and  360-MHz  rays  are  separating  at  a linear  rate,  and  it  is  only  during 
the  periods  when  the  chart  records  show  this  that  the  anadysis  can  be  performed. 

To  ascertain  the  TID  velocity,  three  spaced  differential  phase  sites  could  be 
established.  The  TID  velocity  would  be  cailculated  by  madiing  appropriate  use  of 
the  time  difference  of  the  phase  reversals  at  the  three  sites.  As  stated  earlier, 
these  are  the  times  when  the  amgular  refraction  has  been  reduced  to  zero,  thus 
the  measurements  provide  accurate  information  regarding  the  TID  crest  or 
trough  position,  speed  and  direction. 

Ramge  (r)  to  the  area  of  maximum  density  can  be  measured,  if  phase  data 
as  a function  of  relative  TID  position  among  the  three  sites  is  separated  from 
phase  data  caused  by  ionospheric  refraction.  This  can  be  accomplished  using 
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40-MHi  AMPLITUDE 
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Figure  10.  Changing  Phase  and  Polarimeter  Ramp  Ratio  Caused  bv  Series  of 
TID 


interferometer  analysis  on  the  higher  frequency  (360-MHz)  by  knowing  the  vari- 
ous distances  among  the  three  sites.  However  this  is  the  subject  for  a later 
report,  should  the  expanded  program  materialize. 

Analysis  of  the  data  shows  the  tremendous  influence  a TID  has  over  the 
ionospheric  refraction  of  the  lower  UHF  signals.  In  minutes  it  can  reduce  to 
zero  the  rate  of  charge  of  electron  content  and  temporarily  reverse  the  diurnal 
refractive  process  caused  by  the  sun.  The  reason  for  this  is  the  dominance  of 
the  TIU  gradients  normal  to  the  ray  path  giving  rise  to  large  incident  anglt 
(including  grazing),  as  opposed  to  the  normal,  less  disturbed  vertical  gradients 
with  incident  angles  that  are  much  more  acute.  This  is  in  accordar.ce  with 
Snell's  law,  the  more  oblique  the  angle,  the  lower  the  plasma  frequency,  resul- 
ting in  greater  retraction. 

If  one  accepts  the  foregoing  analysis,  one  must  ask  the  question,  "How  can 

3 

the  ramps  measure  the  electron  density  per  m and,  at  the  same  time,  measure 
the  electron  content  in  a square  meter  column  between  the  earth  and  the  satellite?' 
The  answer  is  that  during  TID  they  cannot,  because  the  formula  for  electron  con- 
tent is  invalid  where  horizontal  gradients  are  present.  This  is  demonstrated  in 
the  following  exercise. 

Shortly  before  the  ATS-6  satellite  disappeared  over  the  horizon,  a 40-MHz 
polarimeter  output  was  fed  on  to  an  adjacent  pen  recorder  channel  for  compari- 
son with  the  differential  phase  data  (Figure  10).  During  this  nine-day  period 
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tiu're  were  sonu'  revolatioiis.  'I'he  ratio  fif  polarimetor  ramps  to 

iliffcrential  phase  ramji.s  varied  during  the  day.  Since  hotli  were  presumed  to  Ije 
a measure  of  the  electron  content  in  a .s.)uare  meter  coIuidji  Ijetween  the  .satellite 
and  the  receivtn-  on  the  (Ground,  except  for  a small  difference  caused  hy  the  pla.sma- 
siJhere,  it  was  expected  that  this  ratio  would  he  nearly  const.'uit.  This  proved  to 
1)0  a false  assumption.  The  two  measured  extremes  are  7 to  1 and  4 to  1 . Shown 
in  the  fifrurc  is  a ratio  of  4.  1:  ramps  to  1 polarimcter  ramp  and  a ratio  of  7 to  1 

at  an  earlier  time. 

In  terms  of  electron  content,  the  formula  for  the  differential  phase  is: 


40  MHz 


where  ' = numljer  of  :iG0'  phase  shifted  cycles  between  the  low  frequency  and 

the  higli  frequency  and  m - ,'160-MIlz/40-MIlz.  The  ramps  displayed  on  tlie  chart 
record  have  been  divided  down  and  represent  12B0'  of  differential  pliase,  making 
one  ramp  equal  to  1.07  ■-  10*  ’ e/ m column. 

The  formula  for  converting  the  polarimcter  measurements  to  electron  content 
is  .\g  = '/(K/f“)  H cos  ft  and  one  ramp  equals  4.7  ■ 10*  ’ e/m^  column,  i'he  dif- 
ference l)etwccn  the  two  measurements  is  important.  The  (.lifferential  phase 
system  measures  the  electron  content  in  the  total  path  (the  ionosphere  plus  the 
plasrnaspherc)  whereas  the  polarimcter  data  is  a measure  of  tlie  electron  content 
onlv  to  the  range  still  influenced  by  the  earth's  magnetic  field.  Thus,  wliere  the 
ramp  ratio  is  4 to  1,  tlie  change  in  electron  content  measures  4.  2li  . 10  e/m" 
column  in  the  differential  phase  system  for  each  3.7  < 10*  ’ e/m*^  column  change 
in  the  polarimcter  system.  This  leaves  a change  of  . 58  >■  10*  ’ o/m^  in  the  m"^ 
column  between  the  magnetic  field  and  the  satellite,  which  seems  to  be  high  hut 
perhaps  reasonable.  However,  when  the  ratio  becomes  7 to  1 , it  represents,  in 

terms  of  electron  content,  7.  4 < 10*'’  e/m'^  to  3.7  y 10*^  e/m*^  which  states  that 

2 

one  half  of  the  electron  change  in  the  m column  is  in  space  beyond  the  iriflu'mcc 
of  the  magnetic  field.  This  is  unrealistic  and  requires  another  explanation  for 
the  inconsistency  in  the  ramp  ratio. 

The  reason  for  the  inconsistency  became  apparent  when  the  TID  records  were 
given  closer  inspection.  Figure  11  is  a 2-hour  graph  of  the  differential  phase, 
polarimcter  ramp  ratio  vs.  time.  Kach  low  point  corresponds  with  the  point  of 
zero  angular  refraction  where  the  phase  reverses  direction  and  the  electron  dif- 
fusion is  greatest  (Figure  10).  The  high  points  coincide  with  periods  of  highest 
phase  velocity  (maximum  angular  refraction).  This  gives  convincing  evidence 
that  the  presence  of  the  TID  is  responsible  for  the  changing  ratio.  It  is  incon- 
ceivable to  think  that  the  TID  could  double  the  influence  of  the  magnetic  field. 


r 
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r lyurc  il.  (hanyiiu;  I’haso  and  I'olarimeter  Hamp  Ratio  Caused  by  TID,  Also 
Sluuviiit;  influence  on  Slab  Thickness 


and  unlikely  that  the  rilJ  would  encounter  sufficient  plasmasphere  electrons  to 
produce  twice  tlie  refraction.  Therefore,  the  conclusion  drawn  from  the  above 
exercise  is  tliut  the  differential  phase  formula  for  calculating  electron  content 
is  not  valid  in  tiie  presence  of  TID,  because  it  produces  horizontal  gradients 
that  dominate  the  physical  interactions;  thus,  the  angular  refraction  is  lietter 
represented  liy  tlie  derived  formula  for  electron  density’. 


').  MMilTIOV  \l.  (. M SKS  OK  IMSCRKPANCIKS  IN  TIIK.  KARADAY/i®  RATIOS 

Sniitti*  ' lias  noted  a diurnal  variation  in  the  ratio  of  Faraday  rotation  T.  to 
electron  content  I of  up  to  10  percent  at  Puerto  Rico.  Kersley  and  Sawbrook'^ 
report  similar  tliscrepancies  in  Great  Britain.  Davies  et  al  are  in  close  agree- 
ment with  .ATS-fi  phase  modulation  and  Faraday  measurements  at  Boulder  showing 
a H"i  percent  diurnal  discrepancy.  They  suggest  that  the  changing  "/l  ratio  .r.ay 
be  due  to  the  diurnal  variation  in  the  effective  height  of  the  ionosphere,  the  per- 
cent reiirosenting  a change  from  580  km  to  1200  km. 

Till-  Hamilton,  Massachusetts  records  were  inspected  for  a changing  diurnal 

'/l  ratio,  but  none  was  found  in  the  six  available  June  records.  However,  it  is 

apparent  in  the  January  12  TKG  profile  comparisons  (Figure  8)  that  in  the  early 

1 1 . Smith,  I).  H.  (Ifl70)  Diurnal  variation  of  the  mean  Faraday  factor  at 
Areceibo,  J.  Geophys.  Res.  7.5:828. 

12.  Kersley,  F.  and  Sawbrook,  I).  J.  (1P71)  Diurnal  changes  of  the  mean 

Faraday  i\1 -factor.  Jfdm  .SatelljUySfudies  Group  Repiirt  1:57.  Institute 
di  Onide  Flettromagnotiche,  Florence,  7taly. 

I J.  Davies,  K.  , Fritz,  R.B.,  Grubb,  R.N.,  and  Jones,  J.  K.  (l!)75)Some 
earlv  results  from  the  ATS-5  Radio  Beacon  Kxperiment,  Radio  Sci. 

10(\o.  9,  10):785. 
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Figure  12.  (Itianging  Flevation 
.\ngle  as  A'l’S-G  Moved  Towanls 
Horizon  Measured  on  IfiO-l'oot 
.-\ntenna 
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Figure  l.'i.  'I'his  Changing  Katio 
Conforms  With  Theoin  ''  = (K/f2) 

Xe  H cos  6 ds  Where  * = Number 
of  Polarinieter  Kamps,  11  = Mag- 
netic I'ield,  6 = Angle  Hetween 
U and  Kay  Path,  9 Was  Increasing, 
Thus  cos  9 Was  Decreasing.  I'he 
Number  of  Differential  Phase  Hamps 
Is  Independent  of  tlie  Magnetic  I 'ield 


morning  tiours  from  0000  to  OGOO  I.T  the  phase  is  increasing  and  decreasing  much 
faster  than  the  Faraday  rotation,  showing  values  that  indicate  al)out  :10  percent  of 
the  electron  content  is  in  the  plasmasphere  at  that  time.  This  would  certainly 
show  a much  lower  ~7l  ratio  and  indicates  a greater  effective  lieight  during  this 
period,  agreeing  with  the  reports  by  Davies  et  al.  ^ 

All  other  varying  ~/l  ratios  were  scattered  through  the  records  and  were 
present  as  reported  above  only  during  periods  of  apparent  horizontal  gradients, 
except  for  a long-term,  changing,  average  ratio  caused  by  the  lowering  of  the 
elevation  angle  as  the  ATS-G  satellite  moved  towards  the  horizon  (Figures  12  and 
11).  Inspection  of  the  TEC  equations  for  both  systems  shows  that  the  changing 
ratio  conforms  witli  theory,  demonstrating  that  only  the  number  of  polarinieter 
ramps  is  related  to  the  magnetic  field  or  more  precisely  to  the  cos  of  the  angle 
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! n'l '.'.f't’t’.  th*'  ■”  ii-  I'ir'ld  :iivl  the  r"i\  I'.'itli.  I'lijs  f-xr-r-rls**  ;il.so  s(‘r",<'S  ii>  I'-i  ' m 

St  r: It  1 I ' I ■ .1 1 : li  t \ I t ! 'I'  I I ■ i :i  1 1 r ipTf. i tiI  ■ Ic.siyti  ;uid  t !u-  l'i(  Id  i • cii  su  n i i - i it  s . 

(I.  I (.11  l\  M l M si.  \l!  nilCKM  SS 

sinci'  tin-  fciri:.  'd'tli<-  ()li;(.sc  data  daritiLt  i'll)  is  consistent  I'/ith  tiit'  fnrn.ula  ■Ic'.t'l- 
oped  fur  etdcuUit n'.L,'  electron  dcnsils  , calculation  of  equivalent  slab  tliickness  is 
also  possilile  if  poUuameter  data  is  uvailalile  to  measure  electron  content. 

Dividmc'  the  electr(.m  content  contained  in  the  Til)  wave  from  the  trouL'h  to 
t!i('  cre-st  (a.s  metesuri'd  bv  tiie  polarimeter),  by  the  electron  density  for  tin'  sami- 
pieriod  (meti.suri'd  b\  the  .1;  method),  \ itdds  the  slab  thickness.  I'he  ratio,  usini.' 
the  parameters  eiven  earlier,  reduces  to: 

(no.  of  i ;i0  polarimeter  ramps)  (d.  h) 

(no.  of  IdliO'  differential  phase  raiiips)  ‘ (sec  z)  ’ 

where  h is  th(>  heittht  and  z is  the  zenith  angle.  It  is  suggested  that  this  ratio  is 
the  sanu'  as  the  ratio  of  'I'KC  to  maximum  electron  densitv  at  the  tinie  of  the 
nieusurt'ment . This  formula  can  be  used  to  assign  values  to  I’igure  1 1,  which 
now  becomes  a very  graphic  portrayal  of  the  slab  thickness  profile  showing  the 
tremendous  itifluence  a TU)  has  on  the  ionosphere.  1 sing  400  km  for  h and  1.  .'i4 

for  sec  z,  it  gives  values  of  slab  thickness  of  120  km  to  202  km.  These  values 

1 4 

are  practicidlv  the  stune  as  the  ranges  reported  by  Tvagi  and  Samayaiulu 

u.sing  the  ratio  'I'KC/max  electron  density.  Thev  reported  a range  of  values 

from  12')  to  200  km  at  Delhi.  However,  vuilues  reported  by  Klohuchar  and  .Allen,* 

with  c.’omhinerl  data  from  two  locations,  are  higher.  Their  ajialysis  on  data  taken 

for  11!  months  in  11'07-1068  yield  a wide  range  of  vedues  with  an  annual  mean 

value  of  201  km.  Thev  do  caution  that  the  values  suffer  somewhat  from  tlie  fact 

that  the  TIIC  data  wore  taken  along  the  oblique  patli  to  tlie  AT.k-:i  satellite  from 

Hamilton,  Ma.ssachusetts,  while  the  N values  were  derived  from  vertical 

max 

incidence  loni.isondcs  from  I'ort  Helvoir,  \Trginia  and  W'allop.s  Island,  Virginia, 

It  aiipe.'irs  that  more  comparison  data  will  he  required  to  iletermine  the 
validity  of  tlie  hypothesis  that  the  change  in  the  TH)  electron  content,  div  ided  b\ 
the  change  in  density,  is  equivalent  to  ' following  exercise  sug- 

gests that  tlieoi-eticallv  tiu'y  are  the  same. 

14.  Tyagi,  I .11.  and  .Samayajulu,  A . V.  (11*00)  .Some  results  of  electron  content 
measurements  from  Karatlav  fading  rif  S-OG  transmissions.  Radio  Sci. 
I:ll2:i. 

10.  KUihurhur,  ■!,  A.  and  .Allen,  R.S.  (1!I70)  .A  first-order  prediction  model  of 
total  electron  content  group  path  delay  for  a midlatitude  ionosphere, 

Al  t HI.  Report  70-1)40  t,  .Air  Force  .Surveys  in  (loophy  sics  No,  222. 
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I'iii-  Klobui'liar  ami  Allen  anals  sis'  ’ showed  that  the  propaiiation  d(dav 

■>  ' • , ■> 

' r In.  1/ f"  \df  where  \ilf  1 I'*'  in  el/ in“  slant  eoluinn. 

Het'ei'i'ine  hai'k  to  the  formula  for  <'leetron  densils  derived  heii'  fur  die  I'll), 

It  is  shown  that  \ ii' f“/  10.  di-  where  \ eleetron  densits  el/m  n-  dr,  tiu.s 

:i-  to.  :;-\/r. 

Sinee  thei-e  is  oidv  oim  dr  representing'  the  hda\  m pro|)apat  i uii  froin  crest 

to  ti’ouch  m th('  Til),  i rrespi-ol  i ve  of  the  cause,  the  two  furmalas  fur  di'  ran  he 

('(|uated,  which  yields  r ' \df/N/m’^.  Thus  it  is  shown  that  such  a relat  lon.-di i p 

liroduees  the  formula  for  the  equivalent  slab  thickness;  however,  in  this  casi’  the 

integration  is  only  durini^  the  I'll)  interval  from  crest  to  trininli,  anil  \ represents 

n 

the  ehanpe  in  ideetron  densitv  for  the  .same  period.  This  is  the  equivalent  m “ 

column  of  h'litrth  r,  with  a uniform  density  of  electrons,  that  is  either  compres.se  1 

■> 

into  or  diffused  avvav  from  tlu'  eiiuivali'ilt  ni“  column  representing  the  total 
electron  content  with  a uniform  ma-vimum  den.sity.  I'lie  i]uestioii  is,  .Are  thev 
reallv  the  same  Iciiftth'.-'  The  analysis  shows  that  at  least  tliiyv  are  m the  same 
ball  park. 


7.  H.M  IHON  COMIM  IN  Tin  1*1. \S\1  V^mil  KK 

( (impannu  I'iuure  10  with  I'iguro  II  shows  that  the  low  points  coincide  with 
times  of  maximum  diffusion  anti  /ero  angular  refraction.  The  indication  is  that 
all  grtidieiit.s  n<.>rmal  to  the  rf  ray  have  momentarily  disappeared,  so  those  are 
the  periods  when  the  ilifferential  phase  formula  for  TMC  is  valid.  Thus,  if 
conqiari sons  are  matle  between  the  I'KC  diurnal  profiles  of  both  systems  to 
measure  the  electron  content  in  the  space  between  the  ionosphere,  innuenced  by 
the  earth’s  magnetic  field  and  the  satellite,  the  Til)  troughs  would  be  accurate 
positions  from  wliicli  to  make  such  measurements.  As  illustrated  in  the  figure', 
the  trough  positions  ptmduce  the  more  stable  and  consistent  values,  indicating 
that  the  averaged  net  TKC  hoyotid  the  magnetic  field  remains  quite  constant, 
even  though  ionospheric  perturbations  are  present  during  this  period. 


H.  \ .><1  NHISK  I’MI  NOMK.NON 

Davies  et  al'  *,  using  l-araday  and  modulation  phase  data  from  an  A TS-fi 
experiment,  have  illustrated  an  interesting  sunrise  phenomenon.  I'hev  have 
shown  that  the  start  of  the  l•■a^aday  TKC  increase  precedes  the  start  of  the 
modulation  phase  TKC  increase,  varying  from  less  than  1 minute  to  1 ,a  minutes 
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It  MinriMi'.  I'hc  n luilula' loii  |il;asr  moasu ronu-nt , Iiki'  the  di ffo rent ial  pliaso 
; ' 1 M si: ri  I : < ■ lit  irl i-il  hiTf,  nu-asiiros  tlic  TIsC  all  tlif  wav  to  tlio  satcllito. 

I 'll-  tiiiii'  ihri'i-i  otu-o  was  all  filuilcil  to  a lowcnn{>  of  the-  effoct ivc  lioiplit  of  the 
ionos[i!ict<  liefoi’e  Ilie  total  content  starts  to  increase. 

I lur  reeor  ls  were  insiieelecl  for  Itiis  sunrise  phenomenon.  There  were  six 
availalile  l erords  with  simultaneous  reeonlinus  of  hoth  (ihaso  and  I 'aradav  . I'fie 
iel.e.  s '.aried  from  a few  .sei'onds  to  two  initiutes. 

< ■' ,i npa ri Hem  of  the  .A  TS-G  differential  pliase  profile  witli  the  A'l'S-d  I 'aradav 
P rofi  le  for  .1  anua  rv  1 J , I >7  . ( I i cure  10)  shows  that  a la  i (>e  port  K>n  rif  the  elect  non 
rnnteiil  wtis  iii  tile  p 1 a s o ; a s ( ill e re  for  several  hours  liefori’  sunrise,  reachiriK  a 
[leak  of  app I'oxni aitelv  ili  (lereent  of  tlie  total  at  Odl'i  I , T . The  data  also  indicates 
a I'airlv  sharp  leni;iri-at ion  between  the  loi.osphere  and  the  pla.sinaspliere  durini.; 
tins  period.  It  IS  reasonable  to  beli<-ve  that  tlie  e.xce.ss  electrons  Were  still 
riowitiL;  into  the  1 rei.oon  at  sutirise,  therebv  lowering  the  effective  heitfht  of  tin? 
ionosphere  tad'ore  new  production  lia  1 actuallv  started,  as  sumtesterl  bv  Davies 
et  al.  ' ' 

In  noth  the  modulation  (ihase  ii o-asureinent s and  the  differential  phase 
; ,easurei;'ent  s,  d,'/  It  ■ 0 occurs  when  the  |)roiiuctton  rate  equals  the  loss 
rate  I,.  ,\t  bloulder,  tin  exae.ole  w;is  ifiveii  sliowmi'  I.  to  be  -'J.  1 fi  v 10^  "r:  ^ser  ' 

at  tlie  I tirtidav  moriniiL:  reversal.  In  accordance  with  the  tibove  statement  at 

, I ■>  -2  - 1 

l,'/dt  0,  2.1h  ■i  10  ■'ill  se<-  lu.st  after  suori.se. 

I'wo  llainiltoii,  Massaidiusett s records  were  tinal\?:ed  for  (he  electron 

produidion  ruti'.  The  analv  sis  \ ielded  values  of  ' 2.10  and  - 2.21  ■ 

12  - - 1 

10  III  "sec  , ayreeinti  tavoralilv  with  the  lioulder,  ( olorado  measurement. 

')  Pii \si  \n,ocnA 

The  lower  channel  in  l n'ure  It)  is  the  analo(>  output  of  a circuit  counting 
the  nundier  of  tO  cvi'les  per  see.  Therefore,  it  is  a measure  of  the  phase 
veloritv  . Note  tfi.it  if  everv  ottn  r peak  is  inverted,  the  trace  approximates  a 
slow  sine  wave  PO  out  of  phase  with  the  trace  of  the  differential  phase  data, 
showini’  It  to  he  the  first  derivative  of  the  phase  analog  (the  summation  of  the 
phase  ramps).  When  it  was  first  recordeil,  it  was  oalihrated  to  re>ud  the  rate  of 
( hange  of  electron  content  pt  r m'  colutnn  between  the  earth  and  satellite.  Ilow- 
e\«'r,  in  view  of  the  new  int»‘rpretation  of  the  phase  data  during  the  periods  of 
Til)  gradients,  it  is  more  accurat<-  to  calibrate  the  data  to  read  rate  of  change 
of  electron  dmisitv  (elertrons  per  m'*)  at  highest  density  region. 
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